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(57) Abstract 

A special mobile radio mobile, vehicle 
or personal location system is described where 
the system operates to locate any one of the 
transmitters in an SMR system. The present 
invention operates by identifying the received 
signals, removing me modulation information 
to reconstruct a demodulated carrier and ana- 
lyzing the carrier signals for multipath distor- 
tion. The multipath disruptions are removed to 
calculate the true incident angle of the trans- 
mitted signal at a base station. A plurality of 
base stations are used to simultaneously cal- 
culate the incident angle of the transmitted 
SMR signal and a central location calculates 
the most probable location of the transmitter 
based upon triangulation. The incident angle 
of the transmitted SMR signal at each base 
station is performed by calculating the phase 
differences of the mcoming voice signal with 
phase ambiguity minimized by using a plural- 
ity of dipole antennas spaced at an irregular 
distance to eliminate phase ambiguity error. 
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nTPF rrnm FINDING AND MOBILE UQCAHON SYSTEM 
FtTR TRIINKFT> MOBILE RADIO SYSTEMS 

5 

PTFT D OF THE INVENTION 
The present patent application relates to radio transmitter 
location in a communication system and in particular to locating a transmitter 
10 to a special mobile radio system 

BACKGROUND OF THE INVENTION 
Mobile tnmking radio is weU known to toe art and is also 
sometimes known as Special Mobile Radio (SMR) systems. A trunked radio 
15 system is a system setup to sei^e a plurality of users to a licensed radio 

frequency band to provide radio or radio-telephone services to base or mobile 
units within a defined geographic market area Typically SMR systems are 
used in the public safety industry for such services as police, fire, and medical 
providers and also in the industrial or land transportation industries, such as, 
20 taxicab services, bus services and the like. Such an SMR system typically 
also includes repeater stations in which the signals from individual users are 
amplified and repeated at various places within the geographic market to 
provide a wider geographic accessibility. 

In the United States, an SMR service provider typically obtains 
25 a license for a band of frequencies within the 800 and 900 MHz frequency 
range from Federal Communication System (FCC). The SMR system is 
described as a trunked system which means that several channels are allocated 
within a given frequency band and the transmitters and receivers use 
automatic switching so that all users can access any one of the allocated 
30 channels which are not currently in use. This form of frequency sharing 

results in a minimum wait or blockage of a message through the system by an 
individual user since the user's transmitter shares all available channels. 

By way of example, the FCC has allocated 250 channels 
located in the frequency band of 809 to 816 MHz. Each channel spacing is 
35 25 KHz and the maximum deviation of a transmit frequency is plus or minus 
5 KHz The required frequency stability is 2.5 ppm for transmitters. The 
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FCC allocation also provides for repeater channels which are located 45 MHz 
above these fiequenries. For purposes of the present discussion, repeater 
frequencies are not required and hence will not be discussed 

Typically the mobile transceivers in an SMR system can be 
5 linked to a public switched telephone network by base station or to cellular 
telqphone systems Access to the public switched network may be from a 
variety of means including DTMF tone signalling for automatic operation 
without human intervention. All the SMR channels are voice channels for 
canying analog voice signals to in-band tone signalling may be used to 
10 interfece to a land-based pubhc switdied telqicme network Links to cellular 
telephone systems has been done m 
Mobile Radio (ESMR) systems. 

Special mobile radio systems to which the present invention 
applies allow the automatic sharing of a plurality of channels by all users of 
15 the system in at blocirii^ Hie system is blocking since 

there are typically more users than there are channels and the channels are 
allocated based on demand Statistically the number of channels allocated to 
a user group is based upon a peak demand usage to minimize blockage of the 
channels due to unavailability of an open channel Of course the greater 
20 number of channels in a special mobile radio system, the fewer times the 
channels are blocked 

Another feature of the special mobile radio system is the use of 
a control channel to allow coherent channel allocation to maximize usage of 
the system The control channel may be a dedicated frequency, a subcarrier 
25 frequency or a sub-audible frequency in selected voice channels in the SMR 
system Allocation and usage of the communication channels is therefore 
negotiated through the control channel in a distributive processing topology in 
which a central control site and the various user sites negotiate the use of the 
channels during handshake transmissions prior to channel usage. Various 
30 types of communications protocol may be enqrioyed to ensure non- 
interference among the various users of the system The control channel 
typically carries control information in a digital format. 
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Prior art Figure 1 shows a typical simplified SMR radio system 
used to communicated with a plurality of cars 101a, 101b, 101c, etc. from a 
control station 103 via cxxtnmunications antenna 102. In this system, all the 
mobile radio transceivers 101a, 101b, 101c share a number of common 

5 frequency channels in ^rtrich an individual t3ransce^(^ may use any one of die 
unused channels at any given time to communicate with the control station 

There is a need in die art, however, to locate the mobile 
transceivers in a special mobile radio system. There is a need to identify and 
locate the mottle transmitters in an SMR system while the mobile transmitters 

10 are transmitting. There is also a need in the art to track the movement of the 
transmitters even if the transmitters change frequency channels over a period 
of time. 

Location of radio transmitters is a technology that has been 
used for many years and in many applications. For radio location, the use of 

15 mobile transmitters and a plurality of fixed base station receivers is generally 
known in the art where the fixed receivers locate the source of the radio 
signal by triangulation of the radio signals. TriangulatiGn is a technique 
where the fixed base station receivers each determine the approximate 
direction of the radio signal from the transmitter. The approximate direction 

20 lines are then drawn on a map where the lines will cross. Since the direction 
lines are approximate, and typically three fixed base station receivers are used, 
the intersecting area between the three lines is usually shaped like a triangle 
since the lines rarely meet at a single point The transmitter is thai located 
within the triangle. 

25 There is a need in the art for a ground-based vehicle or mobile 

transmitter location system which is capable of locating SMR transmitters. 
There is a further need for a ground-based mobile transmitter location system 
whidi is immune from disruption of the radio signals in an urban or irregular 
terrain environment. For these reasons, and for other reasons stated below 

30 which will become apparent to those skilled in the art upon reading and 

understanding the present specification, there is a need in the art for a special 
mobile radio vehicle or personal locator system 
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SUMMARY OF THE INVFNTTfiiNr 
The present invention solves the aboventnentioned problems 
with locking transmitters in special mobile radio systems and other problems 
which will be understood by those skilled in the art upon reading and 
5 understanding the present specification The present invention describes a 
mobile, personal or vehicle location system which locates any one of a large 
number of transmitters in an SMR system. The present invention operates by 
to identify If the location 

of the transmitter requires periodic recalc^^ 
10 any channel frequency changes based on that particular transmitter's use of the 
frequency channels; The received signals are analyzed to calculate the 
incident angle of the transmitted signal at a base station. A plurality of base 
stations are used to simultaneously calculate the incident angle of the 
transmitted SMR signal and a cental station calculates the most probable 
15 location of the transmitter based upon triangulaticHi 

To perform the triangulanon at the central station, the incident 
angle of the transmitted SMR signal is calculated at each base station The 
received SMR signal is detected, the transmitter identity is detennined and the 
modulation information is removed Once the base carrier signal is recovered, 
20 the phase of the received carrier is compared between a plurality of antennas 
m an inegdar spaced array of antennas. The phase differences of the 
incoming carrier signal are compared between the antennas and the phase 
ambiguity is minimized to accurately determine the incident angle of the 
received carrier. 

25 

BRIEF DESCRIPTION OF THE PR A WTNTfy: 
In the drawings, where like numerals refer to like components 
throughout the several views, 

Figure 1 is a prior art description of a typical Special Mobile 
30 Radio (SMR) system. 

Figure 2 is a block diagram of the present invention adapted to 
locate SMR transmitters by tri angulation. 
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Figure 3 is a diagram of the receiving antennas of the base 
stations including a plurality of dipole antennas and a single reference omni 
antenna. 

Figure 4 is a block diagram of a carrier recovery circuit 
5 Figure 5 is a three antenna array for an interferometric direction 

finding system 

Figure 6 is an ambiguity plane plot of the three antenna array 

of Figure 5. 

Figure 7 is the ambiguity plane for the ratios 5:6 of A/2 the 
10 three antama array of Figure 5. 

Figure 8 is diagram showing triangulation of the location of a 
SMR transmitter between three base stations. 

Figure 9 is a diagram of a gnomonic plot of the measured 
bearings of the direction finding algorithm on a sphere. 

15 

HFTATT ED DESCRIPTION OF THE PREFERRE D EMBODIMENT 
In the following detailed description of the preferred 
embodiment, reference is made to the acconpanying drawings which fonn a 
part hereof, and in which is shown by way of illustration specific preferred 

20 embodiments in which the inventions may be practiced These embodiments 
are described in sufficient detail to enable those skilled in the art to practice 
the invention, and it is to be understood that other embodiments may be 
utilized and that structural, logical and electrical changes may be made 
without departing from the spirit and scope of the present inventions. The 

25 following detailed description is, therefore, not to be taken in a limiting sense, 
and the scope of the present inventions is defined only by the appended 
claims. 

The present invention shown diagrammatically in Figure 2 
operates to locate the mobile transceivers in a special mobile radio system 
30 The transceivers may be located in a car 101, in a building (not shown), 

carried by a person (not shown), carried within containers or packages, or any 
number or variety of mobile carriers. In a single transmission, the SMR 
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transmitters are transmitting on a single fiequmcy channd whhin tbe 
channel space. Over an extended period of time, me transnritter may be 
(laiisiiiit a plurality of discrete messages using a variety of channel 
frequencies which may be different for each messages. The location 
5 algorimmsmonnorthecon^ 

frequency changes in me SMR system if a cc«ntinucus tracking of the 
transmitter is required. 

The preferred embodiment of the present invention is designed 
as a separate but parallel operating system used to locate SMR taammtters 
10 wimout any modification^ The present invention 

is designed, in effect, to be an "add-on" system to add the capabilities of SMR f\ 
transmitter location without any modifications to the existing SMR system ^1 
whatsoever. Those skilled in the art will readily appreciate, however, that the 
present invention may be made a partial cff integral part of an SMR system to 
15 save on duplication of hardware and installation costs. 

Typically, the location of the transmitter is determined once but 
for a moving vehicle, periodically updated location irdorrnation may be 
required. Usually, the SMR transmitters are activated only by the operator. 
In a system which is designed to include the present invention as an integral 
part of the overall SMR system, however, the SMR transmitters may be 
polled by a base station or central station to transmit without operator 
intervention so the SMR transmitter location may be determined. This would 
be a "fleet location poll". Also, the present invention may be used as a "panic f\ 
button" or "assistance request" system in which SMR transmitters are used ^ 
25 summon help on operator demand, or automatically summon help such as in 
the case of a crash-activated deceleration switch similar to the deceleration 
switches used to activate vehicle airbag systems. 

As shown in Figures 2, a plurality of base station antennas 
202a, 202b and 202c receive the SMR radio transmissions and the base station 
equipment 203a, 203b and 203c, respectively, each locate the incident angle 
of arrival of the true signal even in the face of selective fading and multipath 
receptions. Several of the most probable incident angles of the received 
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signals are sent from the base stations 203a, 203b and 203c to a central site 
204 for ^angulation. See also Figure 8. 

The location system comprises a plurality of base stations (at 
least two) which transmit signals to the SMR remote mobile units and receive 
5 signals from the SMR remote mobile units. The base stations communicate 
using standard SMR channels carrying voice band information to the SMR 
remote mobile units, or the base stations may simply send alert or 
interrogation commands to activate the SMR remote mobile units. Thus the 
SMR remote mobile units may be activated locally (for example by the 
10 wearer or vehicle) or remotely by me base station. The SMR remote mobile 
units may be used by a person as a communications radio, a radio-telephone, 
a cellular telephone link or as an emergency locator. The SMR remote 
mobile units may also be used to locate the vehicle if stolen, in an accident, 
for vehicle tracking or as an integral communication device. 
15 The SMR remote mobile units communicate over licensed 

channels in an SMR system using a standard SMR infrastructure. SMR 
systems are rnanufactured by vendors such as Motorola and E. F. Johnson In 
the preferred embodiment of the present invention, the identification and 
location determination components of the present invention are designed to be 
20 installed as a parallel system operating completely independent of die existing 
SMR system In this fashion, the present invention may be installed and 
operated independent of the SMR system or it may be made an integral part 
of an existing SMR system 

The transmitted signals from the SMR remote mobile units are 
25 received by the base stations which are then used to locate the SMR remote 
mobile units. An array of irregularly-spaced dipole receiving antennas is 
added to the existing omni receiving antenna at each base station For 
direction finding, the existing omni antenna is used as a reference antenna. 
This is an optional feature since, as will be described more fully below, any 
30 one of the dipole antennas of the array may be used as a reference antenna 
The direction of the received signals from the SMR remote mobile units is 
determined by interferometric direction finding using phase difference 
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comparisons between the various dipoles of the ariay. Unique algorithms 
usmg histogram andysis and spatial fflterm^ 

the received signals at me base stations to eliminate me ambiguity caused by 
signal interference, signal drop out, ambiguily m beating measurements 
received from behind the receiving antennas and multi-path reflections 

The present invention is used in a wide variety of applications 
wherever a special mobile radio system infrastructure already existe The 
present invention may be installed on the radio towers already used to serve 
the SMR system without affecting the operation of such a system. The 
present invention may be used to locate vehicles equipped with SMR devices, 
persons carrying such radios of SMR radios located in homes or bufldings. ^ 
The system may also be used as a "panic button" locator where an operator is 
in need of assistance but may not be able to communicate his identity or 
location. 

* 5 System Overview 

The present invention allows the accurate location of the SMR 
transmitters even in noisy and multi-path environments such as those found in 
urban environments and irregular terrain environments. 

20 Base Station Desip n 

The analysis and decoding of the signals received by the base 
stations from the SMR remote mobile units is done almost entirely in the 
digital domain. The carrier frequencies of the SMR frequencies are down- \ 
converted to a lower frequency in each base station and are then digitally ) 

25 sampled. The digital samples are processed to locate the phase information of 
interest for direction finding. Each base station is constructed with a plurality 
of digital signal processors which enable direction finding of a plurahty of 
simultaneously transrratting SMR remote mobile units transmitting on 
different frequency channels. 

30 The receiving antennas of the base stations are constructed 

using a plurality of dipole antennas 401a, 401b, though 401n (generally 
referred to as 401) arranged in an array, as shown in Figure 3 and a single 
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reference omni antenna 404 The omni antenna 404 is not necessary for the 
preferred emboclimerit of the present invention since any one of the dipole 
antennas of Figure 3 may act as a reference antenna for any of the other 
antennas. The omni reference antenna 404 is used as a reference antenna 
5 because the omni antenna 404 exists as part of the SMR receiving and 
transmitter irnrastructure upon which the present invention relies and is used 
for receiving the SMR voice channel messages from the SMR remote mobile 
units. 

In the prefered embodiment of the present invention, the 
10 direction finding (DF) antennas of each base station uses twenty-four sector 
direction finding dipole antennas 401 arranged m four linear arrays as shown 
diagramrnaticaUy in Figure 3. Each of the four linear arrays consists of six 
dipole antennas and an omni reference antenna The phase between the 
antennas is measured indirectly - each instant of time one of the 24 antennas 
15 is selected and its phase is compared to the omni antenna. 

The switching of the antennas is done every 5 milliseconds 
during the transmission of a particular SMR transmitter. The antennas 401 
are continuously scanned since there could be several SMR remote mobile 
units transmitting simultaneously from an^^ onions. Differentiating 
20 between the various simultaneously transmitting SMR remote mobile units is 
possible. As described above, the simultaneously received signals from 
different SMR transmitters will be on separate channels and the base stations 
are capable of distinguishing and tracking the SMR transmitters 
simultaneously. This is due to the use of a plurality of digital signal 
25 processors operating in parallel to analyze the incoming signals using Fast 
Fourier Transforms (200 Hz bins). The direction finding process, described 
more folly below, is implemented in digital signal processors (DSPs) located 
within each base station and a general purpose computer located in the central 
site. 

30 SMR System Design Considerations 

The present invention must be implemented to operate 
according to the frequencies, modulation techniques and control protocol of 
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the specific SMR system with which the present invention will operate. For 
example, the present ravention wffl 

are using only specific ftequency channels, specific modulation techniques and 
specific identification techniques. For purposes of (lescribing the present 
5 invention, but not by way of limitation, the present invention will be 

described to operate on an SMR system which transmits in the 800-900 MHz 
band using FM modulation of analog \^c« infonnation on the carri^ The 
identity of the transmitter is in an audio subcarrier on the modulated carrier. 
As a first step, me identity of the trarisimtter is easily 
10 determined by eith^^ 

identity (leteimination function of the SMR base station in a separately 
operating component. Hie existing SMR system places the identity of the 
transmitter cm an audio subcarrier although those skilled in the art will readily 
recognize that the ID information could be part of a message preamble, a sub- 
15 audible carrier or transmitted in a parallel control channel. Regardless of the 
technique used, the identity of the transmitter is automatically transmitting and 
automatically determined at the receiver site. 

As a second step, before the location algorithm can begin the 
processing of the phase information of the carrier to determine the incident 
angle of the received carrier, each base station must detect toe presence of the 
carrier so the base stations are synchronized in their operation. Each channel 
in the SMR system will naturally contain background noise and the direction 
finding algorithms should not be allowed to "free run" on noise. Instead, the 
carrier signal must be detected and tracked to signal the start of the direction 
25 finding algorithms in each of the base stations. The detection of the carrier in 
an SMR system is quite simple. In the case of FM modulation of the carrier, 
a simple FM detector circuit is used to detect the presence of an FM carrier. 
The detection of the carrier must be done in each channel where the target 
transmitter is transmitting in the case where the target transmitter is changing 
30 channels. 
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Phase DeteTTTrmalKm of t he SMR Carrier 
Phase conparisons of the received earner are calcdated as 
between a plurality of antennas in an array to determine me incident angle of 
the mcoming carrier signal at each basestation. Since the SMR carrier signal 
5 is severely modulated, phase comparisons cannot easUy be made an the 
modulated carrier. Thus the modulation infbnnatian on Ihe earner must be 
removed before accurate phase comparisons are made. By removing the 
modulation mformation, the carrier signal is treated essentially as a sine wave 
in the calculations described below. 
10 Figure 4 is a block diagram of a carrier recent circuit The 

modulation mformation is removed and the carrier signal is recovered entirely 
in the digital domain. The modulated SMR carrier signal is received at the 
omni antenna and one of the direction finding antennas and the signals are 
downconvertedto a lower hitennediate fiequency suA The 
15 downconverted carriers are digitized by an AID converter at an approximately 
lOOKHz sampling rate. A digital high-pass filter and decimation filters 
further limit the carriers and the conjugate of the DF carrier signal is mixed 
with the omni carrier signal to produce a near-sinusoidal reproduction of the 

unmodulated carrier. 
20 The modulated carrier is also a relatively wide band signal 

which needs to be narrowed The frequency modulated carrier may have a 
frequency deviation of plus or minus 5KHz around a center frequency. The 
phase comparisons between the antennas may not be able to track phase 
changes accurately enough if the slope of the phase changes in the modulated 

25 carrier is too severe. Thus the unmodulated carrier is narrow banded by the 
high pass filters as shown in Figure 4. The direction finding algorithms 
described below are generally immune to some phase distortion and the 
narrow banding of the carrier aids in the accuracy of the overall process. 

The SMR carrier is also subject to selective fading, drop outs 

30 and multipath distortion especially if the transmitter is in a rapidly moving 
vehicle. Thus the carrier may not always be available and the direction 
finding algorithms may be required to operate continuously to update the 
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location of the transmitter as it moves. The steps of detection, identification, 
demodulation and direction finding are performed periodically at preset 
intervals and can also be restarted in the face of ihtenupts 

The determination of the phase of the carrier must first be done 
5 before a comparison of me phase of the carrier can be made between 

antennas. The detection of the phase information of the carrier is performed 
in the preferred embodiment of the present invention by first comparing the 
received phase of the demodulated carrier at the omm reference antenna to the 
received phase of the demodulated carrier at one of the dipole antennas in the 
10 array. A matching filter or autocorrelation function is performed on the two 
received signals until a reasonable match of the phase of me two received 
carrier signals is detenmned This provides the carrier information which is 
then used to synchronize the phase comparisons between the plurality of 
antennas of the direction finding array. 

15 

Antenna Design at the Base Stations 
The DSPs of the base station system compare the exact phase 
arnpHtudes of each of the dipole antennas of the four arrays to the phase 
amplitude at the reference antenna. In the antenna array shown in Figure 3, a 
20 circle surrounding the array is divided into eight sectors. The actual direction 
computation is done in the sectors according to Table 1. 

TABLE 1: Sector Allocation on toe DF Antenna Array 

25 30° to 60° - top and left arrays 

60° to 120° - left array 

120° to 150° - left and bottom array 

150° to 210° - bottom array 

210° to 240° - bottom and right array 
30 240° to 300° - right array 

300° to 330° - right and top array 

330° to 30° - top array 

35 Those skilled in the art will readily recognize that the whole 

direction finding array does not have to be installed in a single location. For 
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instance, each linear array of six dipole antennas could be installed on a 
different fece of roof of a high rise building. In some cases an SMR remote 
mobneunit camot transmii For instance one of me 

sectors feces me sea. In mat case me direction finding array could be 

5 reduced, and only a single fece of 

instaUed at that sector. Usually the interfexmieter does not operate in the case 
where an SMR remote mobile unit is very close to the base station (less than 
1 kilometer). In such a condition, the other base stations -will still perform the 
direction finding and the cenr^ 

10 triangulation 

Far the reasons ecemplified more fuUy below, me spacing of 
the six dipole antennas in each of the four sides of me array is critical to 
maintaining a low degree of ambiguity in direction finding and multi-path 
rejection The base ratios of the base lines was found by simulations with 

15 several relative prime numbers and is selected to be 6 : 10 : 4 : 3 : 5 to 
produce a phase gain PG = 567C The ratio was computed through extensive 
simulations. Although the spacing between the various antenna is critical to 
accurate determination of the incident angle of the carrier, the order of the 
spacing of the antennas in the array is not critical. The order chosen here was 

20 based primarily upon a physical implementation consideration that the array 
frame be mechanically balanced 

Tntfii-ferometric Direction Finding 
The direction finding technique used in the preferred 
25 embodiment of the present invention is performed by measuring the electrical 
phase difference between the antenna of the array caused by the angle of the 
incoming wavefront of the SMR frequency of the channel currently being 
used by an SMR transmitter. In ideal conditions where there is no multi-path 
disruptions of the wavefront, the incoming wave from the SMR remote mobile 
30 unit is corning in angle of 0 from the perpendicular line to the base line on 
which the antennas are mounted 
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Figure 5 is an example of how the interferonaetric direction 
finding technique of the preferred embodimait of the present inv^ 
practiced Figure 5 shows, by way of example and not by limitation, only 
three antennas of the anay. The three antennas are used to deteiinine phase 
differences of the incoming signals between pahs of the antenna As 
described above, the phase difference may also be measured between each 
antenna and a reference antenm such as the omm antenna 

The interferometer direction fmtling (DF) techm 
the directioh of arrival of the wavefront by measuring the phase of the 
mcomingwa^e Win several places in space. This interferometric method 
is a planar ihterferamelric method which assumes that the antenna anay and 
me tr an s mittin g SMR remote mobile unit are all in the same relative plane. " € 

Differences in height of the SMR remote mobile unit and the antenna array 
are ignore and cause little error in the direction finding of the present 
15 invention. 

As seen from Figure 5, the three antenna 501, 502 and 503 are 
irregularly spaced such that distance E>g between antennas 501 and 502 is 
great than distance D l2 between antennas 502 and 503. The distance between 
the antennas is selected to be a multiple of a prime number sequence to 

20 ensure an irregular spacing of all antennas in the array which increases the 
probability mat an electrical phase difference will always be measured 
between some of the antenna members of the array. The distance between the 
antennas is important to minimize ambiguity errors in measuring the phase of > 
me miming wavefront If the antennas are spaced too widely apart, there is ' 

25 increased ambiguity in determining whether the phase difference of the same 
portion of the voice-modulated signal is being compared or whether the phase 
of different portion of the voice-modulated signals are being mistakenly 
compared Since the licensed frequencies of the SMR channels operate in the 
800 MHz band, the wavelength X of the incoming signal is approximately 38 

30 centimeters. 
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Again, referring to Figure 5, the phase difference between the 
two antennas 502 and 503 is know as ^ 
angle of the incoming signal is : 

5 (1) (9^+2^) /2k=Di 2 SIN(B) /X 



whore • 

0 is the angle of the incoming wavefront measured from the 
jq peq)eodicular of the antenna array baseline, 

(Pu is the electrical phase measured between antennas 502 and 
503 

K 12 is an ambiguity factor since the system is capable of 

measuring <j> between ± 71 since the antennas could be 
j5 measuring the phase difference of more than one carrier 

cycle \ 

D 12 is the distance between the antennas 502 and 503, and 
Xis the wavelength of the incoming signal. 
According to this formula, as long as 1} 2 is smaller than A/2, the ambiguity 
20 factor K 12 equals 0. Thus the largest antenna spacing in the array is selected 
to be smaller that 15 altimeters. 

In order to estimate the inherent instrumental error of the 
direction finding interferometer, the above equation can be written in a 
different form: 



25 



30 



(2) 5JJV(6)-(cp 12 +27liC 12 )/( 2 ^i2/^) 

Assuming that all the errors in estimating 6 is due to inaccuracies in 
estimating (p, differentiating the above equation produces: 

(3) o (6) = o (q> 12 ) / (2nD 12 COS{e) /A) 
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where a is the standard deviation 
From this equation, inherent limitations in the system are determined The 
electrical error is divided by the factor 27© 12 fl. which is termed the phase 
gain (PG) which reduces the error by the ratio between the distance and the 
5 wave length Also, the error rapidly grows to infinity as 0 approaches 90°. 
Thus the coverage of the present linear interferometer is hnrited to ± 60° as 
shown by the division of the array into sectors described in above in 
conjunction with Table 1. In order to reduce the error, the phase gain must 
be made as large as possible and each antenna face of the antenna array of the 
present linear interferometer array covers only 1/4 of the space. The antenna 
spacing in the array must have a distance smaller than A/2 in order to resolve 
ambiguities. 

Ambiguities in Interfe rometer Direction Findhi| c r 
15 There are two types of ambiguities in interferometer DF: front 

to back ambiguities (the linear interferometer only covers effectively 1/3 of its 
surrounding area) and multiple possible solutions of the interferometer 
equations, in case where the distance between two antennas is greater than 
V2. This latter ambiguity is caused by measuring the phase difference 
20 between different portions of the voice-modulated signal of the incoming 
wavefront. 

The front to back ambiguity is solved in the present invention 
by using multiple arrays, as shown in Figure 3. In order to cover 360°, the 
present invention uses four sub-arrays of six dipole antennas each, for 

25 covering 90°. The antennas in each array have a vertical beam width of 120°. 
The base station scans all the antennas of the array through the 1:24 antenna 
selector 402 and selects the proper sub-array by comparing the amphtudes 
received from antennas in each array. If two arrays are almost equal in 
amplitude (within 3 dB), then the interferometer process is done on both 

30 arrays and the direction of arrival (DOA) of the signal is computed either 
from both or from the proper array which is selected according to the 
direction of arrival closer to its perpendicular. See Figure 8. 
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The phase ambiguity is solved by using more than two antennas 
in each linear array to compare ihe incoming voice-modulated carrier signals. 
By way of example, but not by limitation, a three antenna array is as shown 
in Figure 5. The distances between the antennas mat are elected as Ey=2*A/2 
5 and D s =3*A/2. Using these ratios, the equation labeled (2) above becomes: 

(4) SIN(O) =<p 12 /27t+jq 2 

10 (5) 5XW(6) =<p 23 /3«+2jq 3 /3 

Figure 6 is an ambiguity plane plot of the three antenna array 
of Figure 5. Figure 6 shows the phase difference (fl 2 plotted against the phase 
difference (Ra as an interferometer ambiguity plane plot Figure 6 shows a 
15 famtty of lines, 601, 602, 603, 604, and 605, parameterized by K, 2 & with 
the same slope created as 0 moves from 0° to ±90°. As derived from 
equations 4 and 5 above, the following formula is used to create the plane 
plot of Figure 6: 

20 (<P2 3 /3*) + 2iC 23 /3= (<p a2 /27t) +JC 12 

Ihe ambiguity resolution algorithm is as follows. The phase 
difference between the antennas is measured as (ft 2 and (fo. On the 
interferometer ambiguity plane of Figure 6, find the line represented by K, 2 

25 and (ambiguity factors) that has the minimal distance, which is the actual 
probability distance, to Cft 2 and (fc. After solving the ambiguity, the present 
invention computes the exact direction of arrival using the full span of the 
array i.e. 2.5X which gives a phase gain (PG) of 57L 

One could select ratios of N:N + 1 of A/2, where N is the 

30 number of antennas in the array. This ratio of N:N + 1 of A/2 may be 

selected to be as large as possible in order to improve the phase gain. But as 
N gets larger, the distance between the ambiguity lines gets smaller and the 
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algorithm becomes vuberable to an error in the proper selection of K, 2 and 
Ka. The result will be a gross error. 

If the measurements distribution of (p is gaussian, then in order 
to compute the gross error probability, the "tail" of the gaussian distribution 
5 from the distance between the ambiguity lines to •'infinity" must be computed 
since the phase is folding aiter ±k This computation described above is only 
an approximation and is almost true for small values of <j((p) which is true at 
high signal to noise ratio (SNR) or high signal to interference ratio. For white 
gaussian noise, a (cp) could be approximated by IN SNR where SNR is the 
10 signal to noise ratio. 

Figure 7 is the ambiguity plane for the ratios 5:6 of A/2. The 
distance between the ambiguity lines became less man half compared to the 
ratio of 2:3 which practically increased the gross error to be impractical. 
Thus, the solution is to add more antennas to the array. Adding more 
15 antennas at different ratios spreads out the ambiguity lines in hyper space 
(multi-dimensional space) with dimension of N-l (N - number of antennas in 
the array) and increases the distance between diem In the preferred 
embodiment of the present invention, it was concluded that in an urban 
environment six antennas are required in a linear array resulting in a five 
20 dimensional plot of Figure 6 (which is incapable of being shown in a 

drawing). The following formula are used to create the plane plot of Figure 7 
which is an example of moving the antenna spacing closer together 



25 



SIN(d) = (A/D 12 ) * ( 9l2 +2nK 12 ) / (2*) = (<p 12 /5it) +2JT 12 /5 



SIN(Q) = (X/D 23 ) * (933+2*^3) / (2it) = ((p 23 /67r) +K 23 /3 
(<p 23 /6it) +K 23 /2= ((p 12 /5ir)+2J!r i2 /5 

Direction Finding Alg orithm 
30 Th e fa" direction finding algorithm relates to measurements of 

a linear antenna array having six dipole antennas per side of a four-sided 
array with a single omni reference antenna (see Figure 3). The calculations 
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arc performed on a six antenna subarray after the sub-array is selected to be 
closest to the SMR remote mobile unit as described above in conjunction with 
Table 1. The full mterferameter formula for a six antennas sub-array is: 



5 (6) 







D 1 




q> 2 






SIN{d) +27t 










fs 




D 5 



















+ 














*5 



10 where 

is the electrical phase measured between the antennas, 
D t is the distances between the antennas 
X is the wavelength 
0 is the direction of arrival (DOA) 
15 is the ambiguity factor, integer value such that, 

-7C<2itDi * SIN(e)/^+2nKi<7i 
e £ is the additive noise (assumed to be gaussian). 
To solve the equation is (6) above, linear array calculations are 
performed to solve the following values: 



20 



WO 95/14936 



PCT/IB94/D0372 



20 



<Pi/ 




2 % 


/2n 






«P 3 / 
/2ti 


id = 


2D /x> 


94/ 
/2lt 




2D /x 


<Ps/ 
/2k 




2D /x 







zj2n 






e 2 /2jt 




}E = 


e 3 /2it 
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(7) 



9 = d X+k + E-V-k=d X+ E 



The results of these calculati ons are as follows: 

the vector elements g> are limited by. | -0.5< <p. < 0.5 j; 
(p is foe measurements vector and d is known; 
-0.5<X<0.5; X is unknown; 

k is unknown; its elements are integer, limited by foe elements 

ofd 

E is foe noise vector that is unknown 
In foe calculations, it is assumed that foe noise is additive such as: 
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1 f° i *j 



15 



£{e,)= 0; i €(1,2,3,4,5) 



2 C 2 
-O 2 

-o 2 



1"J 

i=J-l 
i=j+l 
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2 -10 0 0 
-1-2-10 0 
0-12-10 
0 0-12-1 
0 0 0 -1 2 
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Where : 

<J>, is the electrical phase measured at antenna i. 
(ft = <k - 

25 The solution to this equation is to search among all K e N, the 

specific K that nrimnrizes a weighting function in a wei^ited least square 
procedure (WLS). For every K e N we compute: 
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(9) 

The error will be: 

5 

(10) 

10 Weighting function: 

C(^)= € r ( J K)i?- 1 £(JS) = 
20 

Where A is 5X5 matrix. The rank of A is only 4. In order to 
reduce the amount of real time computation, the following computations are 
performed and the results stored: 

A =A v A (; Ay 

[AvA^eigiA) 
f=a£aJ 

30 Since the rank of A is 4, one of the eigenvalues and thus the 

corresponding eigenvector are 0. Matrices 4X5 and 5X4 are created and 
designated by a ~ over the top of the matrix identifier and thus: 



WO 95/14936 



PCT/1B94/00372 



(13) 

5 



23 



A= F t f 



Using these values in the weighting function produces: 

10 < 14 > tfia-^V*^ 

<P=F<£ 



The full algorithm thus becomes: 
Compute in off-line: 

20 1 . Compute f . 

2 . For ali #e2V compute h= fK> 

3 . Compute the matrix H> 
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Compute in real time: 
1. Compute: 



<P=F<p 



30 

2. For all n e N, confute: 
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cm = 

Find the best three (3) K that mostly minimize C(K). This procedure i, 
done for each hop. 
Compute X(K) for those K selected 

Compute a histogram from all the X(K), using for each DOA its 
distance from the lines (weight) m inverse as weight 
Filter the histogram by passing a boxcar low-pass filter of 3° over the 
histogram in order to cancel small peaks in the neighborhood of a real 
peak. 

Find the highest peak and the second to the highest peak 
Report the results. 

TrianfTiila tion Process 
The input to the location algorithm are reported measured 
directions from several base stations. Each base station reports two directions. 
Each direction is reported with its peak heights which is an indication of its 
validity. The location algorithm projects the direction of the incoming signal 
20 on a gnomonic plane tangent to earth, it filters out the false readings, and it 
computes foe most probable location of foe SMR transmitter. 

Cassini UTM Projection 
The solution of ^angulation on foe earth ellipsoid is a highly 

25 non-linear problem In order to simplify foe problem, we transfer foe 
observed data to a plane called the gnomonic plane as shown in Figure 9. 
The geographic location of the base station on earth is transferred to an X-Y 
position on foe gnomonic plane. Each measured bearing is also transferred to 
the gnomonic plane with a small adjustment. The estimation process is done 

30 on foe plane. The estimated location of the SMR remote mobile unit is not 
transferred back to an earth ellipsoid since the same plane is used for graphic 
display of maps using a Microsoft Windows program or the like. 




4. 
5. 

6. 
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7. 
8. 
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The Cassini UTM projection is widely used for mapping. Its 
advantages are that the mapping are in the Cartesian coordinate system such 
that the E-W and N-S lines are orthogonal. The projection maintains 
approximately both the directions related to the north and the distances at 
5 radius less than 50 M The projection is very convenient for graphic display 
such as in the Microsoft Windows format 

The main distortion of the UIM projection is caused by the 
convergence of the meridian toward me north. There are two 

have to be made: 
10 j Qjn^on of the map north to the 

a bias correction; and 
2. Small correction of the DOA (direction of arrival) of the signal 
according to the location of the base station. 
In actual observations, the correction was actually measured 
15 together with the calibration and measurement of the pointing of the direction 
finding array toward the north The calibration was performed relative to 
known reference point locations in the covered area. To improve system 
accuracy, a reference transmitter having a known exact location may be used 
to correct the system parameters. 

20 

Direction Finding Fixing on a Plane 
The fixing process input are the two DOA's received from all 
the base stations. Those skilled in the art will readily recognize that the 
present invention can operate with more or less DOA's from the base stations. 
25 There are several tests performed on the results received from the base 
stations (BS): 

If the peaks are too low, both measurements are ignored 

If the ratio between peaks is larger than 4, the second peak is 

ignored. 

30 The CS (Central System) computes the fix by using all the combinations from 
the DOA received from the base stations. Note that within one message there 
could be several DOAs. 
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Tie fixing process is an LMS (least mean square) process in 
which the most probable point is found The most probable point is defined 
as the point on the plane that its sum of the squared distances from all the 
DOAs is minimal. The process is as frjllows: 

The distance of a point (X,Y) from a DOA is: 

di lX, Y) ± { Y-YJ cos (0^ - (X-XJ sin (d s ) 

where (X> Y;) are the coordinates of the base station number i. 
dj(X , Y) - distance of (X , Y) from the DOA number i. 
Q - The angle of DOA number i from the norm 
(X , Y) - a point an the plane 
(15) could be written as: 

J^sin^) -3^ cos (Bi) =x sin (6.,) -r cos (0 d ) +d i 

Assuming that M inputs from the base station, (i e {1, >!}), 

equation (16) is written ma matrix form 



10 
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(16) 



(17) 



X ± sinCOJ -Y d cos (6^ 



X H sin(0 M ) -Y M cos(Q M ) 



sin(6 i ) -cosfOj) 



sinfe^) -cos(0„) 



M 



where M - number of bases. As shown in Figure 8, M =3 
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(18) 



10 



15 



X d sintBj) -Y ± cos {Q d ) 



X M sin {Q M ) -Y M cos (0 M ) 
BiniQJ -cos{e ± ) 



sin{6 M ) -cos[B M ) 



20 



Where: 

z is a vector of M elements. 
A is a matrix of M rows, 2 columns, 
c is a vector of M elements. 
D is a vector of M elements. 
The purpose is to find z such that D T D will be minimized Hie 
algorithm is LMS, as described above, and is performed by differentiating g 
relative to c and comparing the result to 0. 
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Com plete Fixing Algorithm 
1 Compute all intersections between all probable DOAs using 
(19). 

2. Perform two dimensional filtering over all intersections, taking 
into account the power of each intersection according to the 
DOAs power that created that intersection The filter size is 
900 X 900 meters. 

3. Find the an intersection from the DOA that are in the filter. 
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Operation Summary 
In summary, and in referring once again to Figure 6, the slope 
of the ambiguity lines is 2/3 and the algorithm finds the closest line to the 
observed point 600. The most probable actual point lies on me line 601. 
5 Figure 7 shows the use of more ambiguijty lines to get better resolution of die 
observed point 701. The actual probable location of the observed point would 
be on the line at position 702 where the slope of each of the ambiguity lines 
is the distance between the antennas. Thus a closer spacing of the antennas 
shown at Figure 7 produces less ambiguity as to the observed versus actual 
10 angle of the direction of arrival (DOA). 

The ambiguity lines shown in Figures 5 and 6 are for three 
antennas. In reality according to the matrices described above, six antennas 
are used in a five dimensional space which is unable to be drawn in a figure. 
In a five dimensional space, the observed point is located, and the nearest 
15 lines in the five dimensional space are located for the most probable actual 
DOA for each based station. By spreading out the spacing of the antennas, 
more ambiguity lines as shown in Figure 7 are used such that antenna spacing 
is best spread out to minimize ambiguity 

In the based station, the observed point 701 representing DOA 
20 0 is compared to the closest three lines in a five dimensional space. A 
histogram is then generated based on the three lines closest to the observed 
point. Then a gaussian distribution of the observed 0 's is computed and 
only the actual points 702 and five dimensional space are picked which He 
within a single standard deviation. A histogram typically produces two peaks 
25 for most probable observed values. Histograms that produce more than two 
peaks are filtered so that only the two strongest peaks are taken as potential 
values. These two most probable actual values of the DOA 0 are then sent to 
the central site for triangulation to determine the position of the transmitter. 
At the central site, the two most probable 0 for the two 
30 histogram peaks around a plus or minus three degree range are recorded for 
two, three or more base stations. Thus 2 n most probable points are received 
by the base station whereas n is equal to the number of stations. Thus if four 
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stations are reporting, two most probable DOA 0, 16 observed points are 
received at the central station A facial Sites: is used with a 900 by 900 
meter space to determine which of the points cluster within the most probable 
location. Within this 900 square meter space, the most probable location is 

5 selected Additional confutation can be made by overlaying the 900 by 900 
square meter space onto a map of an urban area and, for vehicle location, 
determining that the vehicle could only be on a street and not within a 
building. In actual observations, the actual location of the SMR remote 
mobile unit is within a 77 square meter area. 

10 Thus the multipath problem is solved through the use of a 

multidimensional ambiguity space to eliminate observed DOA values which 
are least probable. The use of a plurality of arrays will produce a number of 
observed values within the n-dimensional ambiguity space, some of which are 
due to the actual direction of arrival and some of which are due to multipath 

15 reflections. Since the observed values in the preferred embodiment of the 
present invention are the result of multiple frequencies in multi-channel SMR 
system, observed values and most probable values will be calculated at 
various frequencies if the SMR transmitter uses several channels. Different 
frequencies observed over a period of time may be more or less susceptible to 

20 the multipath problem In addition, for a moving SMR remote mobile unit, 
the doppler shift will be observed but minimized over time since only the 
most probable values due to the strongest signals will survive a histogram 
analysis. Thus the histogram analysis should be performed over a number of 
frequency hops over a longer period of time until the multipath signals 

25 statistically drop out. Thus only the true direction signals will survive the 
histogram analysis and will be passed to the central station for final analysis. 



CONCLUSION 
Although specific embodiments have been illustrated and 
30 described herein, it will be appreciated by those of ordinary skill in the art 
that any arrangement which is calculated to achieve the same purpose may be 
substituted for the specific embodiment shown. This application is intended 
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to cover any adaptations or variations of the present invention. Therefore, it 
is manifestly intended that this invention be limited only by die claims and 
the equivalents thereof 
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WHAT IS CLAIMED IS: 

1 A special mobile radio location system, comprising: 

a special mobile radio having a unique identifier, and having a 
5 transmitter for sending a message on one or more selected frequency channels 
selected from a set of available frequency channels and for coordinating 
selection of the erne or more selected frequency channels on a control channel; 

a first base station having an array of irregularly spaced dipole 
10 antennas operable for receiving the message on the one or more selected 
frequency channels, operable for receiving control information on the control 
channel and capable of determining a first direction of origin of the message; 

a second base station also having an array of irregularly spaced 
dipole antennas operable for receiving the message on the one or more 
15 selected frequency channels, operable for receiving control information on the 
control channel and capable of detemriiring a second direction of origin of the 
message; and 

a central station connected to the first base station and the 
second base station, and operable for 
20 1.) receiving the first direction of origin of the message 

from the first base station; 

2. ) receiving the second direction of origin of the 
message from the second base station; 

3. ) triangulating the first direction with the second 
25 direction to determine the location. 

2. The special mobile radio location system according to claim 1 wherein 
a reference station operates periodically to transmit a reference signal of 
known origin and location to calibrate the first and second base stations and 
30 the central station 
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3. The special mobile radio location system according to claim 1 wherein 
the first base statical and the second base station both include a reference 
antenna such that the phase difference of the one or more selected channel 
frequencies between each dipole of the array of irregularly spaced dipole 
antennas and the reference antenna is determined 

4. The mobile location system according to claim 1 wherein the first base 
station and the second base station both operate to reduce the ambiguity of the 
direction of arrival of the incoming signal by calculating a multidimensional 
ambiguity space having lines indicating the only possible angles of the 
direction of arrival of the signal and detennining the probable points of these 
lines which are closest to the observed value. 



5. The mobile location system according to claim 4 wherein the first base 
15 station and the second base station both generate a histogram of the probable 
values and sending the two most probabl e values to the central station 



10 



6. The mobile location system according to claim 1 wherein the central 
site received the most probable values from the first and second base stations 
20 and applies a spatial filter to only use the most probable values which he 
within a preselected distance of the average location of all most probable 
values, 



7 - A method of locating a trunked mobile transmitter, comprising 

25 the steps of 

a) receiving at a first base station a trunked radio frequency 
signal on a plurality of irregularly spaced antennas; 

b.) removing the modulation information to produce a 
demodulated radio frequency carrier signal; 
30 c.) calculating the observed angle of arrival of the demodulated 

radio frequency carrier signal by comparing the phase difference of the 
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demodulated radio frequency carrier signal between pairs of the plurality of 
irregularly spaced antennas; 

d) detmnining the most probable angle of arrival of the 
demodulated radio frequency earner signal by comparing the observed angle 
5 of arrival to ambiguity lines in an n-dimensianal space where the number of 
anteamas is n+1; 

; e.) repeating steps a - d for each of a plurality of base stations 
to produce a plurality of most probable angles of arrivals; 

£ ) calculating an histogram of the plurality of most probable 
10 angles of arrivals and selecting the strongest most probable values; and 
g.) applying a spatial filter to the strongest most probable 
values from each base station and selecting the location based on the closest 
grouping of the strongest most probable values. 
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